How to model BEC numerically? 
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Abstract 

The new method of numerical modelling of Bose-Einstein correlations ob- 
served in all kinds of multiparticle production processes is proposed. 
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The quantum mechanical in its origin Bose-Einstein correlations (BEC) between 
identical bosons are our basic source of the knowledge on the space-time character- 
istics of the hadronizing objects formed in high-energy collisions [3]. Despite the 
long history of the BEC study the question of their numerical modeling of some 
Monte-Carlo (MC) event generators remains still open. This is because MC event 
generators are (by construction) probabilistic in nature. One is therefore usually 
using some method to enhancing some final configurations of momenta of observed 
secondaries to enhance the observed amout of pairs of identical bosons with small 
mementa differences j3]. Here we would like to advocate different approach using 
observation made j3] that symetrization of the amplitudes of system of identical 
bosons leads to geometrical distribution of the like-particles in the phase-space cell 
they occupy (cf. Table I). There are two possible ways to implement this idea nu- 
merically. The first one (see (21 El for details) is a kind of "afterburner" , which can be 



Table 1: Boltzmann vs BEC 
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used with any MC event generator providing us with distribution of particles. The 
idea is to use the phase-space distribution as given by this MC but to allocate anew 
charges to all particles (conserving the number of +/— and neutrals as given by MC) 
in such a way as to effectively put some particles to the same cell (in which they are 
distributed according to geometrical distribution). Here we would like to go step 
further and present a MC event generator that produces particles according to the 
Bose-Einstein statistics. This can be done in the following way (cf. Fig.l). Particles 



distribution f(E) = exp(—E/T) and allocated to a phase-space cell given by its 
4— momentum p (we call it elementary emitting cell (EEC)). New particles of the 
same charge and energy are then added to this cell according to P = Poexp(—E/T) 
weights until the first failure 1 . The process continues with formation of new EEC's 
until the initial energy W is used up. At the end one corrects for energy-momentum 
and charge conservation. In this way: (i) distribution of EEC's is poissonian and 
(ii) distribution of particles (of the same charge) in a cell is geometrical. This results 
in a Negative Binominal (NB) form of the over all multiplicity distribution. 
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BEC with smearing 
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Fig.l The particles selected according to thermal distribution (upper- left) show no BEC (upper-right). 
Correlations occur only (lower-right) when one allows more particles at the same cell (lower- left). The shape of the 
correlation function C2 (Q) is highly correlated with the spread of the cells in the momentum space. 

It turns out that in order to get proper observed structure of correlation function 
C*2(Q) (and to be able to describe the experimental data, cf. Fig. 2) one has to 
allow for some smearing of momenta of like particles allocated to a given EEC (here 
given, for example, by a gaussian form with width a defining the size of average 
EEC). 

1 It should be stressed that such choice results immediately in Bose-Einstein character of the 
cell occupancy, (n) =P/(1-P) = l/[exp(E/T)/P -l] = l/[exp((E + fi)/T-l] where n = In P Q . 



f(E) = Exp(-E/kT)/E, P = P t exp(-E/kT) 
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Fig. 2 Example of calculations showing that we can fit e + e _ data with some choice of parameters (upper-left) 
and at the same time corresponding multiplicity distribution (upper-right). Lower panels show corresponding 
distributions of the number of EEC's (which has poissonian distribution shape - left panel) and distribution of 
number of particles in EEC (which has geometrical distribution shape - right panel). 



In Fig. 2 we show some example of results for BEC obtained this way (including 
a comparison with experimental data for BEC observed in e + e~ annihilations pQ) 
together with the corresponding results for multiplicity distributions, distributions 
of the number of formed EEC's and distributions of particles in them. Although 
the fit to data is quite reasonably one should treat results presented here only as a 
necessary first proof that such program is in principle possible to be implemented 
and that it is flexible to address (so far only some) experimental data 2 . 

We would like to close with the following remarks. The approach presented here 



2 Our approach is still essentially one-dimensional with transverse momenta entering only via 
their mean value, (jpt), all produced particles are assumed to be 7r( + '~'°) only (no resonances are 
considered either), no corrections for any form of final state interaction is attempted. All these 
factors has to be included before one can use this method to serious description of experimental 
data for all kind of reactions. We plan to pursue such research in the future. 



can be traced to first attempts of describe effect of BEC, only later came space- 
time descriptions which dominate at present and which use symmetrization of the 
corresponding wave functions with the space-time integration over some assumed 
source function (i.e., function describing space-time distribution of points in which 
finally observed secondaries are produced 3 ) jUE]. The shape of this source function 
consists then the main object of investigation. From our approach it is obvious that 
this shape emerges from the necessary spread in the momenta (cf. Fig. 1) and 
therefore is not so much given by the dimension of the source as by the length at 
which particles could still be considered as belonging to a given cell (cf. 0] where 
such distinction has been also introduced but on different grounds). 
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